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Magnetostriction and thermal expansion of the spin-ladder compound piperidinium copper bro- 
mide (CsHi2N)2CuBr4 are analyzed in detail. We hnd perfect agreement between experiments and 
the theory of a two-leg spin ladder Hamiltonian for more than a decade in temperature and in a wide 
range of magnetic fields. Relating the magnetostriction along different crystallographic directions 
to two static spin-spin correlation functions, which we compute with Quantum Monte Carlo, allows 
us to reconstruct the magnetoelastic couplings of (C5Hi2N)2CuBr4. We especially focus on the 
quantum critical behavior near the two critical magnetic fields H c \ and H C 2, which is characterized 
by strong singularities rooted in the low dimensionality of the critical spin-system. Extending our 
discussion in Lorenz et al. [Phys. Rev. Lett., 100, 067208 (2008)], we show explicitly that the 
thermal expansion near the upper critical field H C 2 is quantitatively described by a parameter-free 
theory of one-dimensional, non-relativistic Fermions. We also point out that there exists a singular 
quantum critical correction to the elastic moduli. This correction is proportional to the magnetic 
susceptibility x which diverges as \ ~ 1/ VT at the critical fields and thus leads to a strong softening 
of the crystal. 



I. INTRODUCTION 

By now, the field-induced Bose-Einstein condensation 
of magnonsii^ and its one-dimensional analog^ have 
been observed in many different spin compounds like 
coupled spin-dimer systems jfiifi arrays of coupled spin-1 
chains 9 -*^ or spin laddersi 11 ' 12 Typically, the ground state 
of these systems is a spin singlet, separated by a finite 
energy gap from the lowest triplet excitation. Increas- 
ing the magnetic field, two quantum phase transitions 
are induced. There is a first quantum critical point at 
a field H c i where the spin gap closes, and a second at 
H c2 > H c i, when the fully field-polarized state becomes 
the exact ground state of the system. The phase between 
the two critical fields is gapless and, dependent on the di- 
mensionality, is characterized either by long-range order 
or power-law correlation o 4 ' 13 among magnetic moments. 

Recently, we have reported in Ref. fa high-resolution 
measurements of thermal expansion and magnetostric- 
tion of single crystalline piperidinium copper bromide 
(CsH^N^CuBr.^. This compound has been previously 
identified^ to be a good realization of a two-leg spin lad- 
der. As such, it serves as a perfect model system to test 
experimentally the validity and applicability of theoret- 
ical predictions for quantum critical thermodynamics in 
a controlled fashion. In particular, it has been argued 
that the thermal expansion, a, is an especially useful 



probe in the presence of pressure coupling to the critical 
subsystem because a is then more singular than the spe- 
cific heat. This results in a divergent Griineisen parame- 
ter close to quantum criticality ^ 15 ' 16 ' 17 with an exponent 
characteristic for the universality class of the transition. 
Furthermore, as the thermal expansion is given by the 
entropy derivative with respect to pressure, a oc dS/dp, 
it should exhibit a characteristic sign change near the 
quantum phase transition signaling the accumulation of 
entropy. I- The experimental check of these predictions in 
a controlled model system gives important insights about 
their reliability and applicability in the pursuit of quan- 
tum criticality, especially in materials where the nature 
of the quantum critical point is unknown . 19 ' 20 

In (CsH^N^CuBr^ the presence of two adjacent 
quantum critical points, H c \ and H C 2, indeed leads to 
a rich structure of sign changes of thermal expansion re- 
flecting the positions of entropy extrema in the phase 
diagram. As illustrated in Ref. 13, at low tempera- 
tures a maximum in entropy is located near each of the 
critical fields, H cl / 2 , with an enclosed minimum in be- 
tween. Upon increasing temperature T, the two max- 
ima approach each other and merge at higher temper- 
atures. Moreover, the low dimensionality of the spin- 
ladder results in a diverging thermal expansion^ behav- 
ing as a ~ 1/y/T at the two quantum critical points. 

The compound (CsH^N^CuB^ has a monoclinic 
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crystal structure ((3 ~ 99.3°)2I with the legs of the spin 
ladders oriented along the a axis and the rungs roughly 
(as 20°) along the c* axis of the reciprocal lattice. The 
single crystals used for our study have been grown by 
slow evaporation of a solution of (C5Hi 2 N) 2 CuBr4 in 
ethanol. Typical crystals grow in plates of several mm 3 
with (010) faces and the a and c axes oriented paral- 
lel to the edges^ The measurements have been per- 
formed on a home-built capacitance dilatometer in lon- 
gitudinal magnetic fields up to 17 T for temperatures 
0.3 K < T < 10 K. The absence of three-dimensional 
Need order in (CsHi 2 N) 2 CuBr4 down to temperatures 
Tjy < 100 mK implies a very weak inter-ladder coupling 
giving rise to an extended temperature regime controlled 
by one-dimensional physics.— In this regime, the mag- 
netic subsystem of (C5Hi 2 N) 2 CuBr4 is well-described by 
the two-leg spin ladder Hamiltonian 

N 

T~L = J]] [^J-S^iSj^ + J\\ (SijiSj+iji + Sj^Sj+i^) 

i=l 

- 9 H B H(S? A + Sl 2 )], (1.1) 

where the first and second index specify the rung and 
the chain, respectively. We have identified the exchange 
couplings to be J±/k B = 12.9 K and J\\/k B = 3.6 K, see 
Sec. Ill CI The (/-factor is in fact a tensor and so depends 
on the crystallographic axis along which the magnetic 
field H is applied i^i 

In the following, we shortly review the properties of the 
different ground states of the spin-ladder Hamiltonian 
(jl.ip that appear as a function of magnetic fields At 
zero magnetic field, H = 0, the ground state of the spin 
ladder is a singlet made of short-ranged valence bonds. 25 
Due to the dominant rung interaction Jl/Jii ~ 4, these 
valence bonds can be pictured in a good approximation as 
the singlet states of the spin dimers located on each rung 
of the ladder i 26 ' 27 i 28 The resulting spectrum has a finite 
gap separating the ground state from the lowest-lying 
triplet excitations, and, as a consequence, the spin-spin 
correlations decay exponentially with distance. Upon in- 
creasing the magnetic field, this gap decreases until the 
lowest-lying dispersing triplet excitation touches the sin- 
glet energy giving rise to a quantum phase transition at 
H c \. For higher fields, singlet and triplet excitations hy- 
bridize giving rise to a finite magnetization in the ground 
state. This phase is a Luttinger liquid and has gapless ex- 
citations representing fluctuations of the magnetization 
perpendicular to the applied magnetic field. Finally, for 
even higher magnetic fields the zero temperature magne- 
tization saturates above a critical field, i.e. for H > -ff c2 . 
Here, the excitations on top of the fully polarized ground 
state have again a gap that increases with the distance 
to the transition H — H C 2- 

Whereas in Ref. we focused on the thermal ex- 
pansion and magnetostriction along the c* axis of 
(C5Hi2N)2CuBr4, in the present article we also discuss 
our experimental data measured along the a and b axes. 
Furthermore, we present a detailed comparison with re- 



sults of Quantum Monte Carlo (QMC) simulations of 
the Hamiltonian (|1.1[) for the full magnetic field range 
extending beyond the second critical field i/ c2 . As ex- 
plained in detail in Sec. UH we numerically evaluate for 
this analysis spin-spin correlation functions that are re- 
lated to magnetostriction, as was noted before by Zapf et 
al. in Ref. [2!j. We find excellent quantitative agreement 
between theory and experiment that allows to determine 
the magnetoelastic couplings of (C5Hi 2 N) 2 CuBr 4 and 
their respective uniaxial pressure dependencies. In ad- 
dition, we present a discussion of the quantum critical 
thermal expansion near ff c2 in Sec. IIII Al The effective 
critical theory describing its behavior is known exactly, 
which enables us to directly calculate the thermal expan- 
sion near if c2 without adjustable parameters. Finally, 
in Sec. IIII Bl we note that not only the thermal expan- 
sion but also the elastic moduli obtain a quantum criti- 
cal correction that diverges upon approaching the critical 
fields. This results in a strong softening of the crystal 
that finally triggers a first-order transition in the elas- 
tic system. We present a mean-field discussion of the 
phase diagram with the expected position of the line of 
first-order transitions and the accompanied coexistence 
regions. Finally, we estimate that this effect is too weak 
in (C5Hi 2 N) 2 CuBr4 to be observable because the strong 
one-dimensional signatures will be cut-off by inter-ladder 
interactions before a first-order transition can develop. 

II. MAGNETOELASTIC COUPLING IN 

(C 5 Hi 2 N) 2 CuBr 4 

In the range of temperatures considered here, the main 
contribution to thermodynamics of (C5Hi 2 N) 2 CuBr4 can 
be attributed to the spin system, i.e., the spin ladders. 
As the inter-ladder coupling is sufficiently weak, Tjy < 
100 mK, it can be neglected, and we can approximate 
the magnetic part of the free energy, 

F m = -k B T\ntr{e- A / k * T }, (2.1) 

as arising from an ensemble of decoupled spin-ladders. 
The Hamiltonian thus decomposes into a sum of Hamil- 
tonians, H = ^2 n T~t n , where the index n counts the 
number of equivalent ladders in the transverse plane 
orthogonal to the ladder-leg direction, the a axis in 
(C5Hi 2 N) 2 CuBr4. The Hamilonian 7i n represents a sin- 
gle ladder system as defined in (jl.ip . The free energy 
density is 

where Vb = 859 A 3 is the volume per rung, 22 and N is 
the number of rungs in the ladder. The free energy T m 
deriving from a single spin-ladder Hamiltonian 7i, (jl.ip . 
is given by T m = ~k B T\iYtv{e- n/kBT }. 

The dilatometric properties of (C5H 12 N) 2 CuBr 4 are 
characterized by its elastic energy^ whose low- 
temperature and magnetic-field dependencies we expect 



3 




(a) 



(b) 



FIG. 1: Pictorial representation of the magnetoelastic mode- 
coupling corrections to the elastic energy, (|2.4|l . The wiggled 
line represents strain mj and the dot is the strain-spin interac- 
tion g a . The loop of straight lines in diagram (a) represents a 
two- and in diagram (b) a four-spin correlation function. The 
tadpole diagram (a) results in a force on the lattice respon- 
sible for the magnetostriction and thermal expansion effect. 
Diagram (b) modifies the elastic tensor. 



to be dominated by the spin subsystem. The couplings, 
J a , with a =_L, ||, of the spin-ladder Hamiltonian p.l[) . 
are determined by the exchange integrals computed from 
the electronic wavefunctions for given values of lattice 
parameters. In the presence of strain in the lattice, Uij, 
the lattice parameters slightly change and thus induce a 
variation in the effective couplings, J a — J a (uij). Ex- 
panding this dependence to linear order in the strain, 
J a (uij) ~ J a + g%Uij, with g% = g% = dJ a /dUij\ u=0 , 
we can separate from (jl.ll) a magnetoelastic interaction 
Hamiltonian, 



T^int 

N 

i=l 



(2.3) 



<7™ m S,,i.S.j 



+ .9|"" (S-j,iSi+i,i + 



where summation over the indices n and m is implied. 
The strain u nm fluctuates locally and so depends on the 
site index i. The spin-ladder couples linearly to strain 
and thus acts as a force on the lattice. The lattice re- 
sponds according to Hooke's law and, as a consequence, 
inherits the characteristic temperature and magnetic field 
dependence of the spin-spin correlations in the ladder. 

In order to obtain an effective elastic theory, we con- 
sider the corrections to the elastic energy density, £, 
in second-order perturbation theory in the magnetoelas- 
tic coupling g a . Generally, a correction of order to 
the elastic energy is weighted by a 2n-spin correlation 
function that can be conveniently represented as an n th 
derivative of the free energy F m , Eq. (|2.ip , with respect 
to the couplings J a ; the most important contributions 
are illustrated in Fig. [TJ The resulting effective elastic 
energy density reads 



. 1 1 dF m 

<H M ) = -jUijCijkiUki + y~qJ~ 



la u ij 



+ 0(g 3 ). (2.4) 



The magnetoelastic coupling gives a contribution linear 
in strain iiy represented by the tadpole diagram (a) of 
Fig.ffl The derivative dF m /dJ a acts on the lattice as a H 
and T dependent force that will lead to magnetostriction 
and thermal expansion. In addition, the elastic tensor 
obtains a correction that is second order in the interac- 



tion g, see Fig. [TJd 
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(2.5) 



The tensor c°. fcZ characterizes the elasticity of the crys- 
tal in the absence of spin-strain coupling. Typically, the 
magnetoelastic coupling g a is small and the second-order 
correction in (|2.5[) can be neglected. In the following, 
we will therefore treat the elastic tensor Cjjjy as inde- 
pendent of magnetic field and temperature. As we will 
discuss in Sec. IIIIB| this assumption is very well justi- 
fied in the temperature range of the experimental data 
that we present below. Only very close to the quantum 
critical points, ff c i/ C 2, the tensor Cijki obtains a strong 
dependence on H and T as the second-order derivatives 
d 2 F m /dJ a dJp diverge, resulting in a pronounced soften- 
ing of the crystal. 

Treating the effective elastic theory (|2.4[) on a mean- 
field level, we can minimize the elastic energy (|2.4p . 
d£/duij = 0, to obtain the strain 



1 dF, 



ijkl 



m kl 



V dJ a 



9a = Sal, 



1.) 
a ' 



(2.6) 



where we made use of the inverse of the elastic tensor, 
{c^ r )rjkiCkin m = l(SinSjm + SimSj n ). For later conve- 
nience, we also introduced the dimcnsionless constants 



la 



I ~1\ 1 

-\c )ijki^rg, 



ki 



ijkl 



1 dJ a 



Vd duu 



(2.7) 



u=0 



that quantify the strain dependence of the exchange en- 
ergies; Vd is again the volume per rung. Moreover, we 
defined 



V D dF„ 



I dTr, 



V dJ a N dJ a 



(2.8) 



see Eq. (|2.2p . These quantities can be identified with the 
spin-spin correlation functions along the rung and along 
the leg of the ladder, a =_L, ||, 



S±(T,h) = ^J2(S i , 1 S i , 2 )H (2.9a) 

i 

1 N 

5|, (T, h) = - ^(S^Shu + S,..,S,. 7 , , (2.9b) 



where (6) H = tr{6 e- n / kBT }/tr{e~ n / kBT }. They de- 
pend on temperature, T, and magnetic field, H; the de- 
pendence on the latter only enters via the effective field 
h = gfisH- These correlators will be basic ingredients in 
what follows. They determine the H and T dependence 
of the strain due to the magnetoelastic coupling and lead 
to the magnetostriction and the magnetic contribution 
to the thermal expansion of (CsH^N^CuBr^ 

Uniaxial length changes: The resulting length change 
along some axis n can be obtained by projection of the 
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strain (f2T6l) 



(2.10) 



where n is a unit vector along the n axis. We introduced 
the abbreviation 



1 cV„ 



Vd du 



kl 



(2.11) 



Using (c^^ijki — duki/ddij, with the stress tensor er^-, 
we identify the coefficients 7™ with the derivatives of the 
exchange couplings with respect to uniaxial pressure p n , 



7a 



1 dJ a 
V D daij 



1 dJ a 
V D dp n ' 



(2.12) 



Length measurements on (CsH^N^CuB^ were per- 
formed along the a, b and c* axes. Our detailed com- 
parison between theory and experiment allows us to de- 
termine the coefficients 7" n for n = a, b, c*, which will 
be given below. 

Magnetostriction: The magnetostriction is defined as 
the change of length as a function of magnetic field H 



AL n (T, h n ) 



5Lr. 



(2.13) 



H=0 



The magnetostriction depends on the magnetic field, H, 
only via the effective field h n — g n fj,sH. In the present 
experimental setup the magnetic field has always been 
aligned parallel to the measured length change such that 
the ^-factor also carries the same index n. The g-factors 
in the three directions, the a axis, g a (along the legs 
of the spin ladder), the b axis, gb, and the c* axis, g c * 
(roughly along the ladder rungs), of (CsH^N^CuB^ 
were determined in Ref. to 



g a = 2.06, g b = 2.18, g c < = 2.15. 



(2.14) 



The magnetostriction can be expressed in terms of the 
spin correlators of the ladder, 



AL n {T, h n ) 
L,, 



7 £X>j.(T,& n )+7jfX>||(T,& n ), (2.15) 



where we introduced the abbreviation 

2? ±i |,(T, h) = S ±Al (T, h)-S ul (T,0). 
Similarly, the derivative of the magnetostriction, 



(2.16) 



1 d6L n 



7j 



dS± 
~~dH 



dS\ 



+ 7"^, (2.17) 



can be decomposed into derivatives of the correlators 

Thermal expansion: The uniaxial thermal expansion 
quantifies the length changes of the crystal as a function 
of temperature at constant magnetic field, 



a n (T, h n ) 



1 dSL n 



dT 



11 



H 



d_Sj_ 
~df~ 



dSu 
•7/1 — - 



(2.18) 



It is determined by the derivatives of the correlators <Sj_,|| 
with respect to temperature T . 



A. Experimental results 

Uniaxial magnetostriction of (Cstl^N^CuB^ was 
measured in a range of temperatures from 0.3 K up to 
8 K along the a, b and c* axes, see panel (a), (b) and 
(c), respectively, of Fig. [2l The measurements reveal 
a rich and complex magnetoelastic behavior. From the 
low-temperature data, the critical magnetic fields -ff c i/2 
for each direction can be identified by the kinks in the 
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FIG. 2: (color online). The symbols in each panel show 
the measured magnetostriction AL n /L n of (C5Hi2N)2CuBr4 
along the n axis with n = a,b, and c* for different temper- 
atures. By construction, the magnetostriction vanishes at 
H = 0. The curves at different T are off-set by ±2 ■ 10~ 5 
in panel (a) and (b), and by 1.5 • 10~ 5 in panel (c). The solid 
lines are results from QMC simulations. The insets show the 
fit to the A„ data from which the 7 coefficients in Table|T]have 
been determined, see text in Sec. Ill Dl 
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FIG. 3: (color online). Representative measurements of the 
thermal expansion along the a axis, i.e. the leg direction (left), 
and along the b axis (right) of (C5Hi2N)2CuBr4. The upper 
(lower) panels show data for magnetic fields below (above) 
the upper critical field H C 2- 



30 
20 
L0 

P 

SO 

-20 
-30 
-40 
-50 




a-axis / ( 
c*-axis;\ /*, \ 



b-axis > 
T=330mK 



10 15 
g h H R H/k R (K) 



20 



FIG. 4: (color online). Uniaxial magnetostriction along the 
b axis, Xb see Eq. (|2.17|) . of (C5Hi2N) 2 CuBr4 at different 
temperatures (symbols). The higher temperature curves are 
each offset by additional 1CP 5 /T. The magnetic field is given 
in units of Kelvin, gb^BH /fes, with a g-factor along the b 
axis gb = 2.18. The solid lines are fits using the measure- 
ments along the c* and a axes and relation (|2.19|l with the 
temperature-independent fitting parameters (|2.2U|I . The inset 
compares the magnetostriction along the three directions at 
T ~ 330 mK as a function of the magnetic field H . 



magnetostriction. At small fields, H < H c i, and low 
temperatures, T < Jn/ks, the magnetostriction is expo- 
nentially small as a consequence of the presence of the 
spin-triplet gap. Near H c i the magnetostriction starts 
to grow and saturates near the second critical field H C 2, 
becoming H independent for H > H C 2- For higher tem- 
peratures T > Jy/fes, the kinks near the critical fields 
are smeared out. A saturation plateau at large fields is 
not reached anymore, and the low-field behavior now be- 
comes quadratic, AL n /L n ~ H 2 . The inboxes show the 
derivatives X n , see Eq. pflTf) . for T w 0.58 K < J\\/k B , 
that exhibit clear anomalies at both critical fields H c i 
and H C 2 . As discussed already in Ref . 0, at lowest tem- 
perature the magnetostriction and A c * along the c* axis 
resemble closely the magnetization and susceptibility of 
(C5Hi2N)2CuBr 4f ii respectively. For the measurements 
along the a and b directions, the shape of the magne- 
tostriction derivative is much less symmetric with respect 
to (H c i+H C 2)/2. In particular, A a has only a tiny peak at 
H c i . As will be seen below, this is due to the fact that the 
uniaxial pressure dependence of H c \ almost vanishes for 
a-axis pressure p a due to a partial cancellation of the uni- 
axial pressure dependencies of J± and Jy ; see Eq. (I2.2ip 
and Table |H 

Figure [3] displays some representative data of the ther- 
mal expansion measurements along the a axis (legs) and 
along the b axis. For both directions, we find a rich struc- 
ture in a n with various sign changes, similar to our data 
along the c* axis which have been discussed in detail in 
Ref. [lj- Close to H C 2 there are again clear indications for 
a diverging low-temperature thermal expansion for both, 
a a and u\,. In contrast, however, the features around 
H c i are much less pronounced, which is due the afore- 



mentioned weak pressure dependencies of H c \ for a-axis 
and 6-axis pressure. 



B. Two-parameter scaling 

A remarkable prediction of Eq. (|2.17p is that the mag- 
netostriction curves along the three directions are linearly 
dependent. In particular, the magnetostriction in one of 
the three directions, say, along the b axis can be obtained 
as a linear superposition of the other two, along the c* 
and a axis, i.e. the rung and leg directions, respectively, 

A h (T, h) = Cl A c (T, h) + c 2 A a (T, h), (2.19) 

where c\ and C2 are constants independent of tempera- 
ture and magnetic field. It is important here that one 
accounts for the <?-factor anisotropy (|2.14|) by compar- 
ing the three measurements at the same effective field 
h = g^sH. The sizeable <?-factor anisotropy is reflected 
e.g. in the slight shift of the positions of the high-field 
peaks of X n for the different field directions, as is illus- 
trated in the inset of Fig. [|] 

Fig. [4] confirms that the relation (|2. 19|) is obeyed by 
the experimental data. The symbols in Fig. 2] show the 
measurements along the b axis at different temperatures 
and the solid line is obtained from the magnetostriction 
along the c* and a axes and Eq. (|2.19[) with the fitted 
values 

ci = -0.36, c 2 = -0.85. (2.20) 



This agreement serves as an experimental proof that the 
thermodynamics is fully grasped by a model Hamiltonian 
with only two pressure-dependent coupling constants. 



C. Determination of the couplings J± and Jm 

To perform a quantitative comparison of the dilatomet- 
ric experimental data with theory, we first need to deter- 
mine the two coupling constants, J± and J11, of the spin- 
ladder Hamiltonian (|l.ip . For this purpose, we take the 
critical fields as an experimental input and derive Ji and 
J|| from their second-order strong coupling expressions^ 



9VbH c2 = J± +2 Jm. 




(2.21) 
(2.22) 



The formula for H C 2 is exact in all orders in J\\£& 

To obtain the values of H c \ and H C 2 we make use of 
the temperature scaling of the pronounced peaks near 
the critical fields in the derivative of the magnetostriction 
along the c* axis, A c * (H), see Fig.EK- As the temperature 
is lowered, the peak positions are expected to approach 
the critical magnetic fields. From the known universality 
class of the quantum phase transitions, cf. the discussion 
in Sec. lIIH we expect the magnetic field value of the peak 
positions to scale linearly with temperature close to crit- 
icality, and this is indeed confirmed by the experimental 
data as shown in Fig.JSJ; and d (circles). From the extrap- 
olation we obtain for the critical fields along the c* axis, 
= 6.8 T and H( 2 = 13.9 T. These values correspond 
to critical effective fields, h c i/ 2 = 9c*fJ>BH^, 2 , 



ftd/fc fl = 9.8K, h c2 /k B = 20.1 K. 



(2.23) 



The measurements along a and b yield practically (within 
0.1 K) the same value for h C 2> but due to the small peaks 
around H c i they do not allow for a reliable quantitative 
determination of the lower critical fields. 

Another feature of the peaks in A c * (H) is the steepen- 
ing of its left and right flank near H c \ and H C 2, respec- 
tively, upon lowering T. As a consequence, two A c * (H) 
measurements closest in temperature show a crossing 
point on these flanks whose position approaches the crit- 
ical magnetic fields as the mean of their temperatures 
vanishes. The extrapolation of the crossing point po- 
sition confirms the extracted value of as shown in 
Fig. 0: (stars) . The corresponding test for the H^ 2 value 
is not possible due to the enhanced noise in the experi- 
mental A c * data at higher magnetic fields. 

From the critical fields of the c*-axis measurements 
we obtain the coupling constants by using Eqs. (|2.2ip 
and (|2.22p that are valid up to second order in the small 
parameter J\\/J±. We get 
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FIG. 5: (color online). The derivatives of the magnetostric- 
tion, An of (C5Hi2N)2CuBr4 measured along c* (panel a) and 
along the a and b axis (panel b) for different temperatures. 
Panel (c) shows the linear-T scaling of the positions of the 
maxima of A c * (dots) and of the crossing points (stars) near 
H c i and (d) the positions of the maxima of A n near H C 2 for the 
different axes, n = a,b,c*; see text for further explanation. 



These values are used for the Quantum Monte Carlo sim- 
ulations of the spin-ladder Hamiltonian (|1.1[) presented 
in the following. 



D. 



Quantitative comparison of magnetostriction 
with theory 



J ± /fc s = 12.9K, J||/fc s = 3.6K. 



(2.24) 



With the values for the exchange couplings (|2.24p the 
spin-ladder Hamiltonian describing (CsH^N^CuB^ is 
determined. We computed with the Quantum Monte 
Carlo (QMC) algorithm^ the spin correlation functions 
that enter the dilatometric quantities. In Fig. [6] we 
show the numerically evaluated correlators and T>± 
of Eq. (|2.16p as a function of the effective magnetic field, 
h = g/isH , for a series of different temperatures. 

The rung correlator V± in Fig. [6Jd increases monoton- 
ically with increasing h. This can be simply understood 
in the strong-coupling limit, J± 3> Jii. The rung dimers 
then form an almost perfect singlet at h = such that 
a finite magnetic field augments the triplet component 
to the ground state and thus enhances the ferromagnetic 
correlations. At low temperatures, the rung correlator 
T>± effectively measures the density of triplets and at- 
tains a shape very similar to the magnetization curve of 
the ladder, shown in the inset of Fig. [Ht>- In zero field, 
h = 0, the correlator T>± vanishes by construction while 
at h — > oo it saturates to a value close to one. The devia- 
tion from one is due to the in-chain quantum fluctuations 
at h = that reduce the antiferromagnetic correlations 
along the rung. We can obtain the saturation value from 
the following expansion for the ground state energy per 
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h/k B (K) 

FIG. 6: (color online). Correlators T>± 1 \\, see Eq. (|2. 16|) . of 
the spin ladder Hamiltonian with the parameters (|2.24|l , 
computed with QMC as a function of the effective magnetic 
field h = q^bH (in units of Kelvin) for different tempera- 
tures. The simulations are performed always in the limit of 
large number of sites N ^> Jx/(fc_eT). The inset shows the 
magnetization per spin of the ladder. Note the similarity be- 
tween the rung correlator and the magnetization at low T. 

rung 31 

(2.25) 

For the parameters (|2.24p we get for the saturation value 

1 Fi F 1 

2?x(0,fc)|^ oo = z -^«0.96 (2.26) 

in agreement with our numerical findings in Fig. [6J3. 

At low temperatures, the leg correlator 2?|| of Fig.[6ji is 
instead a non-monotonic function of the magnetic field. 
The existence of a minimum can be understood in terms 
of competing Resonant Valence Bonds42. Below the first 
critical field, the ground state will consist of singlet bonds 
that for J± Ju are mainly formed along the rungs. For 
h > h c i , the ground state acquires a small triplet compo- 
nent, that is reflected in the presence of rare rung-triplet 
states. In the presence of these triplets, stronger singlets 
can be accommodated along the chains enhancing effec- 
tively the antiferromagnetic correlations along the legs 
such that first decreases. However, increasing h even 
further the density of triplets grows, their singlet screen- 
ing clouds start to overlap, and ferromagnetism finally 
prevails leading to a sign change in the function. 

The saturation value at T = can again be derived 
with the help of the ground state energy (|2.25|) , 

^ll(M)L_ = i-g|« 0.75, (2.27) 



axis x 10 5 x 10 5 d In J ± /dp n d\n.J« [ /dp n 

a 4.2 7.7 20%/GPa 133%/GPa 

b -8.2 -6.0 -40%/GPa -104%/GPa 

c* 13.5 -1.1 65%/GPa -19%/GPa 

TABLE I: Magnetoelastic couplings 7^ = {dJ a /dpn)/VD, 
see Eq. (|2.12[) . obtained from a fit of the magnetostriction 
of (C 5 Hi2N) 2 CuBr4 to the QMC data as shown in Fig. [2] 
and the resulting uniaxial pressure, p n , dependencies of the 
exchange couplings J± and Ju . The changes of J± and Jy 
under hydrostatic pressure are given by the respective sums 
of their uniaxial pressure dependencies. 



that coincides with the numerical value in Fig [6k,. 

In Figs. [H we show the comparison between the ex- 
perimentally measured magnetostrictions AL n /L n and 
the Quantum Monte Carlo data. With the help of the 
correlators T>± and T>\\ computed for T = 580 mK as 
a function of H and the relation (|2.15p . we determine 
the mixing coefficients, 7" and Y\\ °f Eq. (|2.12[) . by a fit 
to the experimental data. The resulting values of 7™ 1, 
are listed in Tabled] One can check that the coefficients 
in Table U are consistent with the relation (|2.19|) using 
(|2.20p . Remarkably, these coefficients determined at one 
temperature yield parameter-free predictions for all the 
magnetostriction curves measured at other temperatures. 
As shown in Fig. [2j the agreement between the theory 
and experiment is perfect in a range of T of almost one 
decade. The experimental curves at T < 400 mK have 
not been calculated because our QMC algorithm could 
not produce reliable data at such low temperature. 

The two correlators contribute roughly equally 

to the magnetostriction along the a axis (leg direction) 
and the b axis; the corresponding ratio of mixing coeffi- 
cients is of order one. Along the c* axis (rung direction), 
however, the magnetostriction is dominated by T>± ex- 
cept for a small 8% admixture of the leg correlator. This 
agrees with the naive expectation that squeezing the lad- 
der along the rung mainly affects J± . As a consequence, 
the magnetostriction AL C * /L c * inherits the characteris- 
tic shape of the correlator T>± that resembles at low T 
the magnetization of the ladder as a function of H . This 
also explains the similarity between the susceptibility 11 
of (C5Hi2N)2CuBr4 and the derivative A c *, that shows 
two characteristic peaks located close to the critical fields 
H c x/2- In Ref. 0, we neglected the 8% admixture of the 
leg correlator and interpreted the thermal expansion data 
along the c* direction solely in terms of a pressure depen- 
dence of J±. In this approximation, the thermal expan- 
sion a c * can be related to the H derivative of the entropy, 
a c * oc dS/dH, that allows to understand in simple terms 
the various sign changes of thermal expansion^ 



III. QUANTUM CRITIC ALITY 

Near the two critical fields, H c i and H C 2, the lowest- 
lying triplet excitation or the spin-wave excitation be- 
come gapless, respectively. Both transitions are de- 
scribed by a Hamiltonian of free non-relativistic fermions, 
c(x), representing the corresponding cxcitation a 24 i 33 



dx (a 



h 2 d 2 x 

2m 



+ r)c(x) (3.1) 



where the mass m differs for the two transitions and the 
control parameter, r oc H — H cl / 2 , measures at T = 
the distance to the quantum critical point. The resulting 
critical contribution to the free energy density is 




H=14.5T(exp.) 
critical r—0 




0.25 0.5 1 2 
T(K) 



T(K) 



FIG. 7: (color online). Quantum critical thermal expansion 
near H C 2 of (CsHi2N)2CuBr4 (symbols) along the a and b axis 
and comparison with the parameter-free theory (13. 7[) (lines). 
The critical fields along the three directions as extracted in 
Fig. 01 are given by B.%. = 14.5 T and H b c2 = 13.7 T. 



U(r,T) 



v D h 



''-(k B Tf' 2 T 



k B T 



(3.2) 



with the volume Vd per rung and the distance a between 
the rungs. The scaling function is given by 



T{x) = f ^ log [1 + exp (-y 2 

J0 7T ' 



(3.3) 



A. Critical thermal expansion near H C 2 

In the following, we will focus on the critical behavior 
near the upper critical field, H c2 . The ground state of 
the system for H > H c2 is particularly simple: all spins 
are fully polarized by the magnetic field. This allows to 
determine the parameters exactly^ 



J|,a 2 ' 



r = gfj, B H- (J X + 2J,|) 



(3.4) 



We will consider the predictions for the uniaxial thermal 
expansion a n and compare it to experiment. The uni- 
axial length change close to H C 2 can be obtained from 
(12. 6|) by approximating the magnetic free energy with 
its critical part, F rn /V ~ f CI . Moreover, the dominant 
contribution originates from the sensitivity of f CT upon 
small variations of the control parameter r so that we 
can neglect df cr /dm as it only gives sub-leading correc- 
tions close to the quantum critical point. So we obtain 



(3.5) 



From this expression, we derive the critical uniaxial ther- 
mal expansion close to H C 2, a" = (d6L" /dT) / L n , 




T'{x)-xF'{x) 



k B T 

(3.6) 



where the primes indicate derivatives of the function 
!F{x), (|3.3p . Formula (|3.6p is the limiting behavior of 



thermal expansion near the quantum critical point at 
H C 2 . It is expected to describe asymptotically the exper- 
imental data in the range |r|/fce, T <C J\\/ks = 3.6 K. 
At criticality, r = 0, the thermal expansion diverges with 
temperature as l/v 7 ?,— 



l n \H=H C 



= c 



(71 + 2^)^, (3.7) 



with the prefactor C = (V2 - l) C(l/2)/(2 v / 2^ : ) « 
—0.12066. The next-to-leading order correction to the 
result (|3.7p originates from two-magnon scattering pro- 
cesses that lead to a temperature-independent contri- 
bution to a n at H = H C 2- The strain dependence of 
the mass m in (|3.4p and higher-order derivatives of the 
magnon dispersion, on the other hand, yield corrections 
to (|3.7|) that are suppressed by a relative factor ksT/ J\\. 
All these corrections to scaling are neglected in the fol- 
lowing. 

The mixing coefficients 7"^ entering (|3.6p have been 
already determined in Sec. [n] and are listed in Table U 
This allows us to perform a parameter-free calculation 
of a n near the quantum critical point at H C 2- Fig- El 
shows a comparison of the uniaxial thermal expansion of 
(C5rIi2N)2CuBr4 with formula (|3.7p (lines) on a double- 
logarithmic scale. Note that the critical field H C 2 dif- 
fers for the three crystallographic axes, as it is given by 
H£ 2 — h C 2/ '(SnMs) where h c2 /kB = 20.1 K and the g- 
factors are listed in Eqs. (|2.14p . Measurements of a a were 
performed at the critical field iT" 2 = 14. 5T and are dis- 
played in Fig. Ek,. The data nicely approaches the asymp- 
totic 1/Vt behaviour at lowest temperatures. Panel (b) 
displays «h(T) along the b axis for two fields closest to the 
critical field. The saturation of the curve at H = 13.7 T 
at lowest temperatures suggests that the critical field is 
slightly higher than the value H^ 2 = 13.7 T which follows 
from Eq. (|2 . 23[) with g-factor gb = 2.18. The deviation is, 
however, within the error bar for the critical field value 
of ~ ±0.05 T associated with the analysis in Fig. [5ji. 
As expected, in both panels the corrections to scaling 
behaviour of a n are still sizeable in the considered tem- 
perature range. 



9 



B. First-order transition driven by quantum 
criticality 

The above considerations were based on the assump- 
tion made in Sec. [Til that the elastic tensor (|2.5[) is inde- 
pendent of magnetic field and temperature and, in partic- 
ular, is such to ensure the stability of the crystal. Here we 
show that this assumption breaks down in the immediate 
vicinity of the quantum critical points of the spin ladder. 
Close to the quantum phase transition, the spin ladder 
is described by the Hamiltonian (|3.1[) . Near criticality, 
the most relevant magnetoelastic interaction derives from 
the strain dependence of the control parameter r. As r is 
just a certain function of the exchange couplings J a , we 
can repeat the reasoning of Sec. [II] to derive an effective 
magnetoelastic interaction Hamiltonian, 



dx c' {x)c(x)g l J Uij (x) 



(3.8) 



where the coupling g^P = dr / duij\ u=0 measures the strain 
dependence of the control parameter. In second-order 
perturbation theory in g, this interaction then leads to 
an effective elastic tensor of the form, 



Cijkl 



-ijkl 



(3.9) 



The strong temperature and magnetic field dependence 
enters via the susceptibility Xcr that is obtained from the 
critical free energy 



Xc 



a 2 / c 

d 2 r 



V D h^/k^T 



T" 



(3.10) 



Note that for the spin ladder, r oc H — H cl / 2 , such that 
X cr indeed coincides with the critical magnetic suscepti- 
bility of the system. Moreover, Xcr is positive so that the 
magnetoelastic coupling correction in (|3.9|) reduces the 
elastic moduli. In fact, the susceptibility Xcr diverges - 
like the thermal expansion - as 1/ \/T at criticality, r = 0, 
resulting in a strong softening of the crystal until it be- 
comes unstable sufficiently close but still away from the 
quantum critical point. The lattice is then expected to 
undergo a first-order transition that preempts quantum 
criticality. 

The singular contribution to the elastic moduli in (|3.8p 
is attributed to the low dimensionality d = 1 of the crit- 
ical system, i.e., the spin ladder. Generally, in the pres- 
ence of a strain coupling to the square of the order pa- 
rameter, a divergent contribution to elasticity is expected 
when 



u(d + z) < 2, 



(3.11) 



where v is the correlation length and z is the dynamical 
exponent of the quantum critical point. If this criterion 
is fulfilled the lattice becomes unstable before the quan- 
tum phase transition is reached. The condition (|3.11|) 
is related to the corresponding criterion for compressible 



classical critical systemsi 34 i 35 ' 36 i 37 There, it is known that 
a second-order transition is preempted by an instability 
of the lattice if the specific heat exponent is positive, 
a = 2 — vd > 0. From this, the criterion for compressible 
quantum critical systems (|3.1ip is obtained by replacing 
the spatial dimension d with the effective dimension d+z. 

Compressible critical systems have been studied ex- 
tensively in the past, see Ref. [38j and references therein. 
Renormalization group (RG) treatments of such systems 
yield run-away RG flow that is interpreted as an indicator 
for a first-order transition. For isotropic media, where the 
elastic tensor is characterized by just two moduli, com- 
pression and shear, this conclusion is borne out by ex- 
plicit calculations^ Interestingly, the arising transition 
turns out to be governed by long-range interactions and 
is therefore mean-field like. An analysis of compressible 
critical models for materials with lower crystal symmetry 
is in general, however, rather challenging i 38 i 39 ' 40 

Here, we do not aim for a complete description of 
a possible first-order transition in the monoclinic crys- 
tal structure of (CsH^N^CuB^ triggered by quantum 
critical fluctuations. For an estimate of the location of 
first-order transitions in the phase diagram, we consider 
instead the elastic properties on the level of a simplis- 
tic mean-field theory. Neglecting phonon excitations and 
crystal anisotropics, an effective Landau potential for the 
macroscopic volume change u = AV/V can be derived 

V(«) = ^f-u 2 + f ct {r + gu,T) (3.12) 

where K c g is an effective bulk modulus and with f cr given 
in Eq. (13 . 2|) . The first argument of / cr is the control pa- 
rameter expanded in first order in u with the coupling 
constant g = dr/du\ u=Q . For conditions of fixed hydro- 
static pressure, minimization of the potential V yields the 
value of volume change u m in(fo, T) at given tq and tem- 
perature T. At r = T = 0, the critical part behaves as 
/cr ~ \u\ v ( d+z '&(— u) and dominates over the quadratic 
part of the potential if the criterion (|3.11[) is fulfilled; for 
the spin ladder we have u{d+z) = 3/2. Due to this strong 
non-analyticity of the Landau potential (|3.12|) near the 
apparent quantum critical point, the value of w m in jumps 
in its vicinity and thus avoids the strong critical fluctua- 
tions by prohibiting the argument rg + gu m m to approach 
zero. 

The resulting line of first-order transitions and the co- 
existence region in the phase diagram are shown in Fig. [8] 
The line of first-order transitions terminates in a second- 
order endpoint, which we numerically evaluated to be 
located at 

r*/e = 0.0205 ± 0.0003, k B T*/s = 0.0155 ± 0.0005. 

(3.13) 

4 a 2 2m 

The characteristic energy scale Eq = ,% ^"ti is sup- 

D off 

pressed by the fourth power of the coupling g and 
the temperature T* is therefore expected to be small. 
We expect the second-order transition at the endpoint 
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order transition line (solid) is embedded in a coexistence re- 
gion which is bounded by the dashed lines; the critical fluctu- 
ations at T = are gapped (gapless) on the right (left) hand 
side of the transition line. 



(r*,T*) to be particular: generically, it should be mean- 
field like and as such it is not accompanied by critical 
fluctuations i 37 ' 40 

In order to get the order of magnitude for the transi- 
tion temperature T* we estimate the energy scale So near 
the critical field H C 2 of (CsHijN^CuBr^ Away from 
criticality, the coupling g obeys g = J^r = —K e g-^-r = 
KqsVd'y; with the expression for the mass m, Eq. I]3.4p . 
we then obtain eo = — — P — — ■ Using for the modulus the 

estimate K c g ~ 10 GPa (assuming a rather soft material) 
and a 7 value, 7 ~ 10~ 4 , of the order of the ones in Table|T] 
we get the rough estimate en/fes ~ 10 _5 K. The small nu- 
merical prefactor in (|3.13p suppresses the critical temper- 
ature further by two orders of mag nitude, T* ~ 10~ 7 K. 
This small value of T* suggests that the phenomenon of 
fluctuation-induced first-order transitions close to quan- 
tum criticality is irrelevant for (CsH^N^CuBr^ and can- 
not be observed. In particular, the above discussion ap- 
plies only in the temperature regime where thermody- 
namics is dominated by the one-dimensional physics of 
the spin ladders. At criticality, the inter-ladder coupling 
J3d leads to a dimensional crossover at a temperature 
T^d from Id to 3d quantum critical behavior)^! 41 - which 
weakens the non-analyticity in the effective Landau po- 
tential for the strain u, and Eq. ([3.12)1 ceases to be ap- 
plicable for temperatures T < T%d- Moreover, quantum 
criticality in the 3d regime is not strong enough to have 
the criterion ([3. lip fulfilled. In the framework of the 
3d dilute interacting Bose gas with anisotropic hopping, 2 
the crossover temperature, T$di can be estimated to be 
of the order of the inter-ladder coupling, T^d ~ ^3d/fcs- 
From the measured Neel temperature near (H c i+H c2 )/2, 
Tjv — 80 mK, 22 it follows that is at least three orders 
of magnitude larger than our estimate of T* . This means 
that in (CsH^N^CuB^ the divergent correction to the 
elastic tensor (|3.8p is cut-off by inter-ladder interactions 



before any modulus vanishes. Nevertheless, the strong 
softening of the crystal close to the critical fields should 
give pronounced signatures, e.g., in ultrasound measure- 
ments. 

The derivation of the quantum critical correction to the 
elastic tensor (|3.9|) was based on very general arguments. 
In particular, the prediction of the correction being pro- 
portional to the magnetic susceptibility \ carries over to 
any magnetic-field driven quantum critical point. Such 
a correspondence between \ and the elastic moduli has 
been indeed found by Schmidt et in the spin-dimer 
system NH4CUCI3. This compound exhibits a series of 
quantum phase transitions as a function of magnetic field 
that are ascribed to the successive polarization of nearly 
decoupled spin-dimer subsystems ^ 44 



IV. SUMMARY 

The magnetostriction and thermal expansion of 
(C5Hi2N)2CuBr4 are described to a remarkable level 
of accuracy by the two-leg spin-ladder Hamiltonian of 
Eq. (|l.ip . with the two exchange couplings J± and Jn 
given in Eq. (I2.24p . The presence of only two characteris- 
tic energy scales leads to an interdependence of the mag- 
netostriction measured along different crystallographic 
directions. This predicted dependence is confirmed by 
the experiments, as shown in Fig. [4) and justifies 1 - 1 - the 
disregard of other possible terms in the Hamiltonian, e.g. 
diagonal or ring exchange couplings. As a consequence, 
the contribution of the spin subsystem to the elastic prop- 
erties are captured by two static nearest-neighbor spin- 
spin correlation functions, 2 - which we calculated numer- 
ically with QMC. Their contribution to magnetostriction 
along some arbitrary axis n is weighted by effective mag- 
netoelastic coupling constants that we denoted by 7" n, 
see Eq. (|2.12)l . A detailed fit to the experimental mag- 
netostriction data, see Fig. [51 allowed us to determine 
these coupling coefficients for the three crystallographic 
directions a, b and c* of (CsH^N^CuBrit. Their values 
are given in Table [TJ 

In particular, we find that along the rungs of the lad- 
der, i.e. the c* axis, the magnetostriction is dominated 
by the strain dependence of the rung coupling J±. This 
explains the resemblance of the magnetostriction along 
the c* axis with the magnetization of the ladder.— This 
correspondence was already exploited in Ref. [lj to ex- 
plain the three consecutive sign changes of the thermal 
expansion along c*, that are observed at lowest tempera- 
tures as a function of magnetic field, in terms of entropy 
extremal 

In addition, we analyzed the thermal expansion near 
the upper critical field H C 2- Sufficiently close to the 
quantum critical point at H C 2, the spin-ladder Hamil- 
tonian (jl.ip becomes equivalent to a model of free non- 
relativistic fermions with known microscopic parameters, 
see Sec. IIII Al This enables us to derive a parameter- 
free analytic formula for the thermal expansion along 
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the three crystallographic directions. At criticality, H = 
H c2 , we observe a 1 / VT divergence of thermal expansion 
down to our lowest temperatures. As we discussed in de- 
tail in Ref. Il41 this strong singularity is rooted in the low 
dimensionality of the critical subsystem. 

Finally, we predict a correction to the elastic moduli 
of (C5Hi2N)2CuBr4, whose magnetic field and tempera- 
ture dependence is governed by static four-spin correla- 
tion functions of the spin ladder. Typically, this correc- 
tion is small, and it is indeed negligible for our analysis 
of magnetostriction in the measured temperature range. 
However, close to quantum criticality, the magnetic cor- 
rection to the clastic moduli becomes proportional to the 
magnetic susceptibility x- As the susceptibility diverges 
at criticality as x ~ I/'n/T, this should lead to a strong 
softening of the crystal rendering the elastic system in 
principle unstable at sufficiently low temperatures, see 
Sec. IIII Bl We argued that this phenomenon is in fact 
generic for quantum critical systems if the strain couples 
to the square of the order parameter and, in addition, 
the criterion v(d + z) < 2 is fulfilled. Our estimates 
show, however, that in (CsH^N^CuB^ this fluctuation 
induced first-order transition is not realized as the small 
inter-ladder coupling will cut-off the strong singularities 
before an instability can develop. 



Whereas the strong singularity of x is particular to 
the quantum critical point of the spin ladder, a correc- 
tion to the elastic tensor proportional to the magnetic 
susceptibility is a generic feature of magnetic-field driven 
quantum criticality. Such corrections have been observed 
in ultrasound experiments on NH4CUCI3 by Schmidt et 
alM 

As an outlook, we remark that the precise knowledge 
of the magnetoelastic couplings of (CsH^N^CuB^ as 
given in Table U will be especially useful in the quest 
for a theoretical explanation of its thermal transport 
properties^ 
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